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Tj" ! Abstract 

(N- 

A PT-symmetric optically active medium that lases at the threshold gain also acts as a 
complete perfect absorber at the laser wavelength. This is because spectral singularities of 
PT-symmetric complex potentials are always accompanied by their time-reversal dual. We 
investigate the significance of PT-symmetry for the appearance of these self-dual spectral 
singularities. In particular, using a realistic optical system we show that self-dual spectral 
q-( singularities can emerge also for non-PT-symmetric configurations. This signifies the existence 

of non- "PT-symmetric coherent perfect absorbing (CPA) lasers. 

<N; 

The recent interest in the role of PT-symmetry in optics was initiated with the theoretical works 
\q '. predicting unusual properties of PT-symmetric optical lattices pQ and the experimental works [2] 
^ ! studying the occurrence and consequences of exceptional points corresponding to the spontaneous 
^ \ breakdown of PT-symmetry [3]. Exceptional points have been a subject of both theoretical [I] and 
O ! experimental [5] research for more than two decades, but these early studies did not consider systems 
^ i possessing PT-symmetry. The study of optical realizations of PT-symmetry and its breakdown 
^ | that is carried out in [TJ |2] focuses on the behavior of electromagnetic waves propagating in a 
waveguide consisting of balanced gain and loss regions such that the gain/loss properties of the 
^ \ system varies along a normal direction to the longitudinal axis of the guide as depicted in Figure [TJ 
Another recent development is the discovery of an optical realization of the mathematical concept 
of "spectral singularity" [6] in the form of zero- width resonances [7]. The physical model used in 
[7J to characterize these resonances also involves a PT-symmetric optical waveguide (Figure [2]), 
but unlike in [H [2] here the gain/loss properties of the optically active region of the guide changes 
along its longitudinal axis. This particular waveguide configuration was initially proposed in [8] as a 
physical model possessing PT-symmetry. In [9] it was used to address the problem of constructing 
a metric operator [10] for a pseudo- Hermit ian Hamiltonian [H] with a purely continuous spectrum. 
This in turn paved the way for the discovery of the optical realizations of a spectral singularity 
[TJ, as the latter appeared as a singularity of the metric operators associated with a delta-function 
potential with a complex coupling constant [T2] . 

Ref. [13] shows that the zero-width optical resonances, that we refer to as optical spectral 
singularities, can also appear in non-PT-symmetric waveguides. Ref. [H] uses an optical system 



s E-mail address: amostafazadeh@ku.edu.tr, Phone: +90 212 338 1462, Fax: +90 212 338 1559 



1 



Figure 1: (Color online) Schematic view of an active dielectric waveguide with gain/loss properties 
changing along a normal direction to the longitudinal axis of the guide. 

Guide axis 




Figure 2: (Color online) Schematic view of a waveguide containing an optically active region whose 
gain/loss properties change along the longitudinal axis of the guide. 



consisting of a single gain region to show that the mathematical condition for the occurrence of 
an optical spectral singularity (OSS) coincides with the physical threshold condition for lasing. 
Therefore, OSSs correspond to lasing at the threshold gain. 

At the values of the parameters of the system that a spectral singularity emerges, the transmis- 
sion coefficient |T| 2 as well as both the right and left reflection coefficients, |i? r | 2 and \R l \ 2 , diverge, 
while the wave number k (and consequently the wavelength A := 2ir/k) remain real [TJ. This is 
because the (complex) transmission and reflection amplitudes, T and R r '\ are related to the entries 
Mjj of the transfer matrix M of the system according to 

M 22 ' M 22 ' M 22 ' 1 ; 

and a spectral singularity appears whenever we can satisfy 

M 22 = 0, (2) 

using a real value of A, [7J. In this case the system amplifies the background noise and emits 
radiation of wavelength A from both ends of the scattering region. The time-reversal of this phe- 
nomenon corresponds to an optically active region absorbing incident coherent radiation of identical 
amplitude and phase from both its left- and right-hand sides. This happens at wavelengths where 

m 

M u = 0, (3) 

and provides the basic principle governing the behavior of a coherent perfect absorber (CPA), also 
known as an antilaser 
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In order to make ([3]) more explicit, we invert ([T]) and use the fact that the transfer matrix has 
a unit determinant HE] to obtain 



M 



R l R r R r 

T — 

T T 



R l 1 
Y T 



(4) 



In light of this equation, we can express (JBJ) in the form 

T 2 = R l R r . (5) 
This corresponds to a zero of (one of the eigenvalues of) the S-matrix of the system, 



S := 



T R r 
R l T 



while (EJ), that is equivalent to 



Bf,Er,T-> oo, 



(6) 



(7) 



marks a pole of (one of the eigenvalues of) S, [7J. 

It is easy to see that equations (|3j) and fljTJ) are not contradictory. This shows that there may exist 
a wavelength A at which one of the eigenvalues of the S"-matrix diverges while the other vanishes. 
This corresponds to the situation that a spectral singularity accompanies its time-reversed dual. 
We refer to such a spectral singularity as a "self-dual spectral singularity." An optical system 
supporting a self-dual spectral singularity with wavelength A emits electromagnetic radiation at 
this wavelength unless it is subjected to incident radiation of the same wavelength and identical 
amplitude and phase from both sides. In the latter case, the system acts as a CPA. Therefore, it 
is called a "CPA-laser" [17]. 

Self-duality of spectral singularities is a characteristic feature of one-dimensional PT-symmetric 
scattering potentials. To see this we recall that the transfer matrix of such a potential satisfies, for 
real values of k, M* = M -1 , [To]. Therefore, 



Re(M 12 ) = Re(M 21 ) = 0, 



M, 



22- 



In particular, any zero of M 22 is also a zero of M u . This shows that spectral singularities of VT- 
symmetric potentials are self-dual. Furthermore, one can easily show, using fll]) and (jSJ), that for a 
PT-symmetric scattering potential, 



e 2i(ipi-ip t ) _ e 2i(ip r -ip t ) _ _^ 

\T\ 2 ±\R l R r \ = 1, 



(9) 
(10) 



where {pi,(p r ,{p t are respectively the argument of R l , R r ,T, i.e., e 



i-pi 



T/\T\, and the unspecified sign in (TTOj) stands for e 



i(ifl+ip r -2ipt) 



symmetric complex generalization of the standard continuity relation, 
recently reported in [T9j . 



= R l /\R% e^ r := R r /\R r \, 
Equation (JTU]) is the VT- 
\T\ 2 + \R r > l \ 2 = 1, that is 
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Figure 3: (Color online) Schematic view of an infinite planar slab of optically active material that 
consists of two layers of equal thickness and different loss/gain properties. 



In this article, we use a simple yet sufficiently rich optical system that allows for a clear assess- 
ment of the importance of PT-symmetry for realizing a CPA-laser. Figure |3] shows a schematic 
description of this system. It consists of a two-layer infinite planar slab of optically active mate- 
rial. The layers have equal thickness, L/2, and different complex refractive indices, ni and 1x2- We 
examine the linearly polarized time-harmonic electromagnetic waves propagating along the normal 
direction to the slab. If we use a Cartesian coordinate system whose positive x- and z-axes are 
respectively along the polarization and propagation directions, the wave equation satisfied by the 
electric field E(z,t) = e~ tult ^(z)e x reduces to the Helmholtz equation 



*"(z) + fc 2 n(z) 2 *(z) 



0. 



'ir, 



where u and k := uj/c are respectively the angular frequency and wave number, e x is the unit vector 
pointing along the positive x-axis, ^> : IR — > C is a continuously differentiable function, and 



rii for < z < 0, 



n z 



n 2 for 
1 for 



< z < k 



(12) 



z > 



Our aim is to examine the consequences of imposing ([2]) and 0. 

In [201 121] we derive a formula for the transfer matrix of complex potentials vanishing outside 
a closed interval. In order to make use of this formula, we introduce the dimensionless quantities 



x :- 



z 1 

I + r 



&:=Lk, i/j(x) :=*(Lx-f), v(x) := & 2 [l - n(Lx - f ) 2 ] 
and express (TTT]) as the time-independent Schrodinger equation: 

- ip"(x) + v(x)xjj = tfijjix), 

where 



2)1 



£ 2 (1 - n 2 ) for x e [0, \ 
v ( x ) = { £2(! _ n |) for x e 

for x i [0,1]. 



(13) 



(14) 



(15) 
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Note that v is a PT-symmetric potential provided that n* = n 2 . 

It is not difficult to show that the transfer matrix of the potential ( fl~5i) has the form [21~] : 



M 



2i& 



where 



- e -^(r 1+ - 2r 2H 
e «( ri _ - 2r 2 _) 



i>'(l)±z 



(16) 



and 0^ are the solutions of (I14p in the interval [0, 1] that satisfy the initial conditions: 
<f>'i(0) = 02(0) = 1, and 0' 2 (O) = 0, [20]. They have the following explicit form. 



(17) 
>i(0) = 1, 



<f>j(x) :-- 



A. e iM!X + B . e -iSlnix for x e 1^ 
C . e i&n 2 x + D . e -i^ 2 x for x e ^ 1 j j 



where 



A 1 

£>i 

C 2 



ni 



:= - 1 + _ 
ni 



-4, 







l\ . 


f 1 — — ) cos ai - 




sin ai 




\n 2 


my 



n 2 



1 H I cos ai — i I — H ] sin ai 



n 2 ni 



cos ai + 



mi sin ai 



(19) 
(20) 
(21) 
(22) 

-3 - "2"- (23) 

Notice that the change of variable: z— > x — z/L + 1/2 involves a shift in the original of the z 
axis that changes the form of the transfer matrix of the system. It is easy to compute its effect and 
show that the transfer matrix associated with the coordinate system centered at z = is given by 



n 2 



D 



2 ■= 



COS di 



ztii sin Oi 
n 2 



and 



M 



2ift 



r^(T 1+ - 2r 2H 

rv - sr 2 „ 



Ti+ 

-e^rv 



(24) 



According to this equation, spectral singularities and their time-reversal dual are respectively given 
by the real values of & satisfying 



rv 
r 1+ 



o, 

2r 2+ 



o. 



(25) 
(26) 



In particular the system supports a self-dual spectral singularity and acts as a CPA-laser provided 
that we can satisfy both these equations for a real A. Substituting (j!8p - (122 j) in ( II 7p . we find after 
a lengthy set of calculations the following simple expressions for Ti_ and Ti + — 2r 2+ . 



ri. 

Ti. 



2nin 2 

2r 2+ = 



n + n+ sin a+ + n_n_ sin a_ + i(n 2 + cos o + — v?_ cos a_) 



2nin 2 



n+n + sin a + + n_n_ sin a_ 



i{n + cos o H 



n cos a_ 



(27) 
(28) 
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where 



n± := tii ±n 2 , 



ft± := nin 2 ± 1 = -(n 2 - n 2 ) ± 1, 



a± := Oi ± a 2 



n±ft 



In light of (E7D and (EHD, we can easily reduce ([25]) and ([26]) to 



n + n + sin a + 



-n_n_ sin cl 



n, cos a + = n_ cos a_ 



(29) 



(30) 



Squaring both sides of these equations and using the the trigonometric identity, sin 2 a±+cos 2 a± = 1, 
and the last equation in f[29]) . we obtain 



cos 



cos 









2 J 












r) 







where 



s :-- 




n_n_) 2 — (n + n+) 2 



n_n+, 



n+n- 



(31) 
(32) 

(33) 



Eqs. ( 13T1) and ( 1321) are a pair of complex transcendental equations involving two complex and one 
real variables, n± and respectively 

First, we examine the special case that the refractive indices rii and n 2 have the same real 
part, say n Q , i.e., they correspond to different amounts of gain/loss achieved using the same active 
medium. Denoting the imaginary part of tv,- by Kj, so that tv,- = uq + iKj, we find n_ = — k 2 ) 
and n + = 2n + i{k\ + k 2 ). Inserting these relations in ( 13TT) gives 



[k\ — K2) cosh 



±[2n + i(ki + k 2 )]s. 



(34) 



Because the left-hand side of this equation is real, so must be its right-hand side. However, notice 
that for non-exotic material 1 < uq < 5, and \kj\ are at least three orders of magnitude smaller 
than n , i.e., \kj\ ^ 10 _3 n . Therefore, we can approximate s by expanding it in power series in kj 
and ignoring the quadratic and higher order terms in the latter. This gives 



n 2 n + l 



2m (Ki + k 2 ) 



n 



Now, we insert this relation in (I34p and demand that its right-hand side be real. The result is 



[ki + k 2 ) 1 + 



0. 



(35) 



Because n > 1, 



1 + 



> 1 



> 1. 



Therefore, f[35]) implies that K\ + k 2 w 0. This means that within the limits set by the physical 
considerations, the condition that the system supports a self-dual optical spectral singularity implies 
that rii ~ n^, i.e., the system possesses PT-symmetry. However, there are more general situations 
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where rii and ri2 have different real parts and we can still satisfy (l3Tj) and (132]) . These correspond 
to non-PT-symmetric self-dual spectral singularities that we wish to characterize. 

We can solve Eqs. (13"Tj) and (I3"2"j) for R and use the identity, cos^ 1 ^) = ±iln(j ± y/$ 2 — 1), to 
express the result in the form: 

2i 

R = ± — \ In |u_| + i [arg(±zu_) — 27rm_] k (36) 



2/ 



[in 


u_ 


[in 


Kl 



(37) 



Here we have introduced 



'n 2 s 2 . n_s /nis 2 



+ 1± , u + := W^- + l±-^, (3* 



ni n+ \/ ni n_ 



and employed the identity In 3 = In I3I + i[arg(3) — 2irm], arg(^) stands for the principal argument of 
the complex number 3 (that takes values between —it and tt), and m and m± are arbitrary integers. 

Notice that although & is a real variable, ( 136]) and (137]) are actually complex equations. In 
particular, ^ is equal to the real part of the right-hand side of these equations, and their imaginary 
part must vanish. Therefore, if we denote by r]j and Kj the real and imaginary part of Uj, for 
j — 1, 2, so that rij = rjj ± inj, and let 

r/± := Re(n ± ) = Re(ni) ± Re(n 2 ) = r]i ±r/ 2 , (39) 

n± := Im(n ± ) = Im(ni) ± Im(n 2 ) = «i ± k 2 , (40) 

V?_ := arg(±iu_), (p + := arg(±z u+), (41) 

we find 

£ = ±2|n + |~ 2 [K+ln|u_| +r/ + (27rm_ -</?_)], (42) 

£ = ±2|n_|~ 2 [K_ln|u+| +77_(27rm + -¥>+)], (43) 

?7 + ln|u_| - hi + (2irm^ - </?_) = 0, (44) 

?7_ In |u+| — K_(27rm + — <^ + ) = 0. (45) 

We can solve for In \u±\ in the last two equations and substitute the result in the first two to obtain: 

ft=±— — , R = ±^ 46 

V+ V- 

Because for realistic situations, 1 < r] + < 10, < |r/_| < 10, \<p±\ < ir and 8. > 10 2 , ( 146]) implies 
that ±m_ ^> 1 and ±sgn(?7_)m + ^> 1, where sgn(r^_) stands for the sign of 7?_. In view of these 
relations, we can redefine m± so that m± > and express ( 1461) as 

^2 (ftrm. TP-) (47) 
^= 2(27rm+TV?+) . (48) 
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In terms of the new m±, (144)) and (1451) take the form 



|?7„| In |u_i_| , , 

27rm + =F ip + 
, ry+ln |u_| 

k+ = ±- . (50) 

27rm„ =F 

Because m ± 3> 1 > |<£>±/(27r)|, we can neglect terms proportional to <p±/(27rm±) in f)4T|) - (jSTfl) . 
This gives 

47rm_ . . 

£ « , (51) 

V+ 

Si ~ ^t, (52) 
|?7-|m|u + l 

«- ~ ±— ^ , (53) 

. ?7 + ln |u_| , r 
k+ « ± '+ 54 
27rm_ 

Equations (15 ip and ( 152]) imply the curious relation: 

|?7_| m + 



77+ m_ 



(55) 



In particular, because \f)-\ < rj + , we have m + < m_. 

Next, we note that for realistic models \k±\ J$ 10~ 3 and r] + > 1. Therefore, if we consider the 
cases where |?7_| > 0.1, we have <C \f]±\- Neglecting terms of order \k±/t)±\ in (138]) . we then 
find 



(T/g + lK^ + l) / [(r7 + + ^) 2 + 4][(r7 + -^) 2 + 4] 
W - 1)(% 2 - 1) V [(»?+ + V-) 2 ~ MV + ~ V-) 2 - 4] 



,£J -¥ + l±^, |u + |^/^± + l±^±, 



These relations together with (15 ip and (|52p allow us to write (153"]) and (I54p in the form 




a 2 rj\ 


M± 




v 2 - 


V- 


cr 2 ?] 2 . 


M± 




2 

'/+ 


V+ 



(56) 



^'Hft+^j- (57) 

where the pair of unspecified signs appearing in the right-hand side of these equations are unrelated. 

The approximate equations (15T|) . (152"]) . (I5E|) and (I57P can be used to express m± and k± in 
terms of 771, 772 and ^ up to sixteen sign ambiguities that we eliminate by enforcing ( 130]) and 
using physical considerations. The resulting approximate values can be the basis for more accurate 
numerical solutions of fl30|) . Table [1] gives the result of three sample calculations. The property 
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Approximate Values 


Exact (Numerical) Values 




3.600 ± 1.180 x 10~ 3 z 


3.603 ± 1.178 x 10~ 3 z 


n 2 


1.500 =F 2.241 x 10" 3 2 


1.498 =f 2.243 x 10" 3 i 


m 


3.600 ±2.524 x 10~H 


3.600 ±2.520 x 10~ 3 z 


n 2 


3.000^2.698 x 10" 3 i 


2.997^2.695 x 10" 3 2 


ni 


3.000 ± 1.372 x 10~ 3 z 


3.000 ± 1.370 x 10- 3 z 


n 2 


1.400 =F 2.429 x 10~ 3 i 


1.398 =F 2.431 x 10" 3 i 



Table 1: Results of some typical values for the parameters of the system that yield a non-VT- 
symmetric CPA-laser operating for 8. = 4007T. This corresponds to L/X = 200. The exact values 
have been obtained by a numerical treatment of ( !30|) after fixing the value of 

that these configurations come in complex-conjugate pairs is a manifestation of the fact that the 
time-reversal transformation maps each solution of the equations determining self-dual spectral 
singularities, i.e., M\\ = M 22 = 0, to a solution of these equations with the same value of the wave 
number As we see from Table [H our simple approximation scheme produces values which are 
in a very good agreement with the exact (numerical) results. 

In conclusion, in this article we have performed an analytic treatment of the problem of finding 
self-dual spectral singularities of an experimentally accessible model. Self-duality is a characteristic 
feature of the spectral singularities of complex "PT-symmetric scattering potentials. However, as 
our explicit calculations show, one can also achieve a self-dual spectral singularity for a non-VT- 
symmetric potential. In terms of the optical realizations of spectral singularities, this signifies the 
existence of a large class of non-"PT-symmetric CPA-lasers. These may be of practical interest, for 
these devices are not bound by the strict requirement of having balanced gain and loss componentsJ§ 

Acknowledgments: This work has been supported by the Scientific and Technological Research 
Council of Turkey (TUBITAK) in the framework of the project no: 110T611, and by the Turkish 
Academy of Sciences (TUBA). I wish to thank Aref Mostafazadeh for his help in preparing the 
figures. 

Appendix 

In order to give an idea of the amount of freedom we gain by relaxing the condition of "PT-symmetry 
in devising a CPA-laser, consider a system that can be modeled using a complex scattering potential 
involving n independent complex coupling constants, 31,32, ■ ■ • , In- We can write it in the form: 
v n{x) = Y%LiUfe~( x ) + *Z)f=m+i3</i"( a; )) where // : R R and ff:R->R are respectively 
even and odd real-valued functions that form a linearly-independent subset of the vector space of 
all real- valued functions defined on R. 

1 This is true for every scattering potential. 

2 A scattering potential involving n independent complex coupling constants has n — 1 degrees of freedom, if we 
demand that it is 'PT-symmetric and has a self-dual spectral singularity. This number is 2n — 3, if we relax the 
requirement of 'PT-symmetry. For a proof, see the appendix. 
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Requiring that v n be PT-symmetric implies that all the coupling constants are real. In this 
case, Mn = M| 2 , and the condition that v n has a spectral singularity, which is necessarily self-dual, 
means to choose these n real coupling constants and the wave number k such that M 22 = 0. Because 
this is a complex equation, it restricts two of the these n + 1 real parameters. Therefore, imposing 
PT-symmetry and the condition of the existence of a self-dual spectral singularity amounts to 
fixing the n + 2 of the initial In + 1 real parameters of the problem, and leaves us with n — 1 degrees 
of freedom. 

Now, we relax the condition that v n is PT-symmetric and only demand that it has a self-dual 
spectral singularity. This means that the coupling constants 3^ and the wave number k satisfy 
Mn = and M22 = 0. These are a pair of generically independent complex equations. They 
restrict 4 of the initial 2n + 1 real physical parameters: Re (3^), Im(3^) and k. Therefore, in this 
case we are left with 2n — 3 degrees of freedom. This is larger than that of the PT-symmetric case, 
if n > 2. For the simple model that we consider in this paper, n = 2 and the number of degrees 
of freedom left after imposing the condition of the existence of a self-dual spectral singularity is 
the same for both the PT-symmetric and non-PT-symmetric configurations. However, even in 
this case the range of values of the physical parameters supporting a self-dual spectral singularity 
can be altered almost at will for non-PT-symmetric configurations. Specifically, we can freely 
choose rji, 77 2 , and & in our approximation scheme and find exact values of rji and 772 in a close 
vicinity of our initial choice for these quantities. For the PT-symmetric case, we can only choose 
77 = 771 and .ft. Therefore, although relaxing the requirement of PT-symmetry does not increase 
the degrees of freedom of the system, it provides an extra freedom in the choice of the range of 
values of one of the physical parameters. For n > 2, non-PT-symmetric configurations (optical 
systems) supporting a spectral singularity (functioning as a CPA-laser) have a larger number of 
free parameters. Therefore, they should be practically more favorable over their PT-symmetric 
counterparts. 
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